from scipy import *
import numpy
n=10
def t0(p):
    return numpy.cos(numpy.pi*(2*p+1)/(2*n))
xx = arange(n,dtype=longfloat)
print xx
#print t0(xx)
print numpy.cos(numpy.pi*(2*xx+1)/(2*n))